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The translational motion of molecular ions can be effectively cooled sympathetically to trans- 
lational temperatures below 100 mK in ion traps through Coulomb interactions with laser-cooled 
atomic ions. The ro-vibrational degrees of freedom, however, are expected to be largely unaffected 
during translational cooling. We have previously proposed schemes for cooling of the internal de- 
grees of freedom of such translationally cold but internally hot heteronuclear diatomic ions in the 
simplest case of ^E electronic ground state molecules. Here we present a significant simplification 
of these schemes and make a generalization to the most frequently encountered electronic ground 
states of heteronuclear molecular ions: ^E , ^E , ^E and ^11 . The schemes are relying on one or 
two laser driven transitions with the possible inclusion of a tailored incoherent far infrared radiation 
field. 
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I. INTRODUCTION 

The cooling and manipulation of neutral molecules has 
become the subject of intense studies in recent years 
and impressive advances have been made. Experiments 
include the successful production of molecular Bose- 
Einstein condensates [ll Hi 12] j the deceleration and trap- 
ping of polar molecules in inhomogeneous fields |l,|a,lail3 
and loading a trap with paramagnetic molecules cooled 
by a He buffer gas 8, 9,]. For the NH radical the 
presence of an unusually large Franck-Condon factor of- 
fers prospects for direct Doppler cooling of a trapped 
molecule |Ty|- 

Molecular ions constitute another class of molecules 
that are very interesting to cool and manipulate. Di- 
atomic molecular ions are, e.g., important constituents 
of interstellar media [ill Il2 | , comets and cool stellar at- 
mospheres including that of the sun |12. \W{ , and there 
have recently been proposals to utilize cold molecules to 
implement a quantum computer [ij, |l5j . 

The cooling of molecules is in general more complicated 
than that of atoms since the ro-vibrational substructure 
of the electronic molecular energy levels normally makes 
it impossible to find a closed optical pumping scheme to 
be used for conventional laser cooling. Molecular ions, 
however, may be very effectively cooled sympathetically 
by loading them into a trap with laser cooled atomic ions 
HE HI EE 111 El- The Coulomb interaction between the 
charged particles provides efficient momentum transfer 
from the initially hot molecular ions to the cooled atomic 
ions. Dissipative cooling of the translational motion is 
hence obtained for both species although only the atomic 
ions are subject to laser cooling. 

One might expect that the ro-vibrational degrees of 
freedom of a diatomic molecule placed in the vicinity of a 
cooled atomic ion would couple to the translational mo- 
tion of the atomic ion, resulting in strong sympathetic 



cooling of these degrees of freedom. In a typical ion trap, 
however, the excitation energy of the translational atomic 
motion in the trap (vibrations in the harmonic trap po- 
tential) is of the order of 1 MHz, which is much smaller 
than typical energies of ro-vibrational excitations (of the 
order 10^^ — 10^^ Hz). The large difference between these 
numbers prohibits that the internal ro-vibrational states 
couple effectively to the external motion of the ions in 
the trap. In the following we therefore assume that the 
internal degrees of freedom relax to equilibrium with the 
black-body radiation (BBR) present in the trap. This 
will happen on a timescale of tens of seconds, which is 
significantly faster than the inelastic collision time in the 
trap which, from Langevin theory, is estimated to be hun- 
dreds of seconds ^^ . 

In Refs. |2ll l2a | we proposed schemes for cooling of 
the rotational degree of freedom of such molecular ions in 
the case of heteronuclear molecules with a ^E electronic 
ground state. The schemes are based on two direct in- 
frared (IR) transitions between the lowest vibrational 
states in the molecule or two Raman transitions cou- 
pling the vibrational levels via a near-resonant excited 
electronic state. In addition to the pumping by the 
external light sources the cooling schemes are assisted 
by rotational redistribution mediated by the BBR. The 
timescale of the cooling schemes are on the order of ~ 60 
s which is shorter than the estimated inelastic collision 
rate with background gas. 

Though most molecules appearing in nature have 
a ^S electronic ground state it is necessary to con- 
sider other electronic states for molecules produced in 
the laboratory, including molecular ions. The by far 
most frequently encountered electronic ground states of 
such molecules and molecular ions are, apart from the 
^S state, the ^H , ^E and ^E states. This includes the 
lighter diatomic hydrides, e.g., FH+ (^H ), BH+ (^E ) 
and 0H+ ('^E ). Such ionized hydrides are attractive can- 



didates for our cooling schemes as they have low reduced 
masses and hence high rotational transition frequencies 
leading to fast rotational relaxation rates which is benefi- 
cial for the timescale of the cooling scheme. By extending 
the schemes to ^H , ^E and ■^E states we have then cov- 
ered all the lighter ionic hydrides and the vast majority 
of other molecules amenable for cooling. We show that 
it is possible to cool such electronic states, though at the 
cost of introducing more laser frequencies in some cases. 
For most molecules, however, the present cooling schemes 
rely on only a single IR laser, possibly assisted by broad- 
band radiation from a far-infrared (FIR) emitter which 
is filtered to optimize the cooling efficiency. 

The present paper is organized as follows. In 
Sec. m we present cooling schemes for the ^S electronic 
ground states. In Sec. IIIII we discuss a model of the 
cooling schemes and present numerical simulations for 
MgH+ (X^E). In Sec. lIVI we present cooling schemes ap- 
plicable to the ^E , •^E and ^11 electronic ground states 
together with numerical simulations of each of the cool- 
ing schemes. A summary of the results is given in Sec.lVl 
In Appendix ^ we have collected the Einstein coeffi- 
cients for the considered molecules and transitions, and 
in Appendix IbI we describe the Honl-London factors of 
interest. 



II. COOLING SCHEMES FOR ^E STATES 

The suggested schemes for ^E states are sketched in 
Fig. n The driven transitions are either Raman transi- 
tions via an excited electronic state or transitions directly 
between vibrational levels. Fig.^a) represents the cool- 
ing scheme of Ref. [2lj in which two Raman transitions 
make a closed cycle through pumping of population from 
the "pump states", {ly = 0,N = 1) and (iy = 0, TV = 2), 
to the excited states, {ly — l,N = 1) and {ly = 1, N — 0), 
respectively, followed by subsequent spontaneous emis- 
sion bringing the populations back to the "pump states" 
or to the ro-vibrational ground state. Here z/ and N 
denotes the vibrational and rotational level respectively. 
Population initially in higher-lying states is fed to the 
pump states through BBR-induced rotational transitions 
within the vibrational ground state. 

It would be advantageous for practical implementa- 
tion to use only a single Raman (Fig. Qlb)) or a sin- 
gle direct (Fig. Q^c)) transition, at the expense of not 
emptying the {v = 0,N — 1) state. Without apply- 
ing other means to limit the pile-up of population in the 
{ly = Q,N = 1) state the cooling efficiency, measured as 
the percentage of population in the ground state, will de- 
crease. One can, however, take advantage of the higher 
frequency of the {ly = 0,N = 1) ^ (jy = 0, iV = 2) 
rotational transition compared to the undesired {v — 
0, N = 0) — > {ly = 0, N — 1) heating transition in- 
evitably driven by the BBR, and apply an incoherent 
source and a high-frequency pass filter to reduce the ra- 
diation resonant with the heating transition while still 



addressing the {p = 0,N ^ 1) ^ (ly ^ O^N = 2) 
transition. Thereby one can obtain the desired deple- 
tion of the {ly — 0,N = l)-population by means of 
incoherent radiation only. As the rate of depletion us- 
ing realistic incoherent sources will be slower than if the 
state was addressed by a laser, it is necessary to design 
the cooling scheme such that spontaneous decays to the 
{ly = 0, N = 1) state from states which are participat- 
ing in the pumping cycle are avoided. This can be done 
by addressing the (i/ = 0, TV = 2) ^ (i^ = 2,iV = 0) 
transition with a resonant, dipole allowed (AA^ — 0, ±2) 
Raman pulse as depicted in Fig. ^b). The pumping to 
the {ly = 2,N = 0) state is then followed by spontaneous 
decays through {ly — 1, N ~ I) to {ly = 0,N = 0) and 
{ly — 0,N = 2) in accordance with the dipole selection 
rules {{AN,Aiy) — ±1) for single photon decays. 

It is shown in Ref. l21l that the Raman-transitions 
in the MgH+ test case T(| are saturated by a ~ 100 
kW/cm^, 10 ns pulse, which is a modest intensity for 
present day laser systems. 

One of our schemes using only a single direct laser- 
induced transition subject to the dipole selection rule is 
shown in Fig.QJc) [23|. The laser pumps the (z^ = 0, A^ = 
2) ^^ {ly = 1,N = 1) transition while subsequent spon- 
taneous decay brings the population back to the pump 
state or to the ro-vibrational ground state. A filtered in- 
coherent source is applied in order to bring population 
from the (z^ = 0, A^ = 1) state to the "pump state". The 
advantage of this direct scheme is that it does not depend 
on the existence of an excited electronic state that can be 
addressed with laser light and also that it requires only a 
single laser frequency. ArII+ is an example of a molecule 
without excited electronic states [2lL |23J | . 

From a practical point of view, a pulsed laser system 
is desirable for the direct scheme of Fig. ^c). The IR 
light could, for example, be generated by difference fre- 
quency mixing of the primary beam of a frequency dou- 
bled Nd:YAG laser and a dye laser pumped with the 
same beam. In the MgH+ case, the wavelength of the 
{v = 0,A'^ = 2) ^^ {i> ~ 1, N = 1) pumping transition 
is ~ 5.9 /im I^J and the Einstein A-coefficient is ~ 20 
s^^. To ensure saturation of the laser driven transition 
we require that the population in the states involved un- 
dergoes at least 10 Bloch oscillation during a laser pulse 
and that the amplitude of each oscillation exceeds 0.9. If 
we assume a detuning of 1 GHz and a pulse duration 10 
ns we find that an intensity of ~ 500 W/mm^ is needed 
to fulfill both requirements. This corresponds to a pulse 
energy of 5 /iJ. Typical nonlinear crystals should be able 
to deliver an energy of ^ 10 fi3 per laser pulse at the 
wavelength required. 

The added incoherent field from a lamp will increase 
the rate of rotational transitions needed for cooling, 
but at the expense of heating the population distribu- 
tion. The spectral distribution of the incoherent field can 
therefore be shaped to maximize the cooling efficiency as 
described in Ref. 1231. 



III. NUMERICAL SIMULATIONS FOR 

^S STATES 



o 
II 
> 



A L 



(a) Cooling scheme of Ref. |M 
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In this section, we present our model of the coohng 
scheme, show the results of numerical simulations and 
discuss the optimal radiation distribution of the incoher- 
ent field. 



A. Rate equations for the population dynamics 

The population dynamics is well-described by rate 
equations giving the change in population of a given state 
via Einstein coefficients and frequency-specific radiation 
intensities. The equation of motion for the molecular 
population Pi in state i takes the form 
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(2) 
represents the populations in vector form with Nmax cho- 
sen so the population in higher-lying rotational states 
is negligible during the cooling process. Aij and -By- 
are the Einstein coefficients describing spontaneous and 
stimulated transitions from energy level i to j. Wiuiij) 
is the cycle averaged radiative energy density present in 
the trap at the resonant transition frequency lo = u>ij, 
between level i and j. In Eq. (^), the first term corre- 
sponds to spontaneous decay from state i to states with 
lower energy, while the second term describes sponta- 
neous decay from levels with higher energy into state i. 
Stimulated emission from the ith state and stimulated 
absorption from lower-lying states is then described by 
the third and fourth term, and finally, the last two terms 
represent transitions due to absorption of radiation from 
the ith state and stimulated emission from higher-lying 
states into the ith state. 

The system of Eqs. (Q is conveniently expressed by 
the matrix equation 



(c) Direct scheme 



FIG. 1: Ro-vibrational states of interest in the cooling 
schemes for ^S states. The cooling concept involves tran- 
sitions between ro-vibrational states driven by Raman pulses 
(solid lines, double arrows in (a) and (b)) via an excited elec- 
tronic state or direct laser pumping (solid line, single arrow 
in (c)) and subsequent spontaneous decays (dashed lines). 
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where K is an {M -f 1) x (M + 1) coupling matrix. 
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B. Calculation of molecular properties 

As seen from Eq. (QJ, it is necessary to know the Ein- 
stein coefficients to simulate the population dynamics. 
For many molecules, the Einstein coefficients are avail- 
able in the literature. If not, they are evaluated numer- 
ically as follows. We use the well-known quantum me- 
chanical expressions for the Einstein coefficients between 
an upper state y„ and a lower state ^^ that are both 
non-degenerate |25| 
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where lo denotes the transition frequency and D the tran- 
sition dipole moment between the states. 



D\::'„, = / ^>iM''^''^„,dT, 



(5) 



with dr denoting the volume element corresponding to 
integration over the complete set of coordinates for all 
particles involved and 
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the dipole operator. 

The equations refer to a laboratory-fixed coordinate 
system so the molecular wave-functions include the rota- 
tional terms. The summation indices in Eq. 10, k and 
Z, refer to the electrons and the involved nuclei, while Zi 
denotes the nuclear charge. 

For degenerate states Eq. Q) is modified to 
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where dipole matrix elements connecting ro-vibrational 
levels, 'Dn,rm are derived in Appendix IBl 
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FIG. 2: Top: Born-Oppenheimer electronic potential energy 
curves of MeH^ (X^E) calculated by Gaussian in a 6-31H — hG 
basis set |31| using Hartree-Fock (HF) theory, M0ller-Plesset 
nth order perturbation theory (MPn), coupled cluster theories 
with singles and double excitations (CCD), and single, double 
and triple excitations (CCSDT). See Refs. IH IH IM HI HI 
for descriptions of the methods. The curves for the MP4, 
CCD and CCSDT calculations are in good agreement close to 
the equilibrium position, 1.65 A, indicating that these meth- 
ods give an accurate description of the problem. Bottom: 
Corresponding dipole moment functions, D^°^{R) of Eg. lBSII 
pointed along the internuclear axis, of MgH+ calculated with 
Gaussian. The MPn and coupled cluster theories largely agree 
around the equilibrium distance, although not as well as for 
the potential curve due to the dependence on electronic wave 
functions rather than eigenenergies. The classical turning 
points for the vibrational ground state are marked on the 
common abscissa at 1.5 and 1.8 A. The result of the MP2 
calculation cannot be discriminated from the MP4 result at 
the internuclear distances of interest. 
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The Honl-London factors, Sj^j^, are tabulated in the 
literature 24, 26, 27, 28] and may be evaluated by the ex- 
pressions given in Appendix B. Both the potential energy 
curve for the molecule and the electronic dipole moment 
function £)™°' (i?) is evaluated with Gaussian 29] . From 
the potential energy curve the ro-vibrational eigenfunc- 
tions, /i,^, are readily found using the Numerov method, 
and the one-dimensional integral of Eq. (jSJ can be evalu- 
ated. We use the Level 7.5 program 30] to perform these 
tasks and to evaluate Eq. {Tj) , leaving us with the desired 
Einstein coefficients. 



1. Einstein coeffictents for MgH^ 

Since translational cold samples of MgH+ have been 
produced in a trap loaded with laser cooled Mg+ atomic 
ions [ig, this molecular ion is the first choice for an im- 
plementation of the presented cooling schemes. To our 
knowledge only a few Einstein coefficients for transitions 
within the electronic ground state have been published 
(2lj]. We have re-calculated the coefficients using the ap- 
proach of the previous section. The potential curves ob- 
tained from Gaussian |29| using various theoretical ap- 
proaches on a 6-311-I--I-G basis set jsj is given in Fig. [21 



together with the corresponding dipole moment functions 
in the molecular center of mass system. The potential 
curves show convincing convergence and our derived vi- 
brational transition frequencies and equilibrium distance 
agree with published data within 1.5% 24] . To compute 
the accurate electronic dipole moment function is more 
challenging, since this requires accurate electronic wave 
functions. Generally the MUer-Plesset fourth-order per- 
turbation theory [33 and the coupled-cluster theories |3g 
are reliable for the task. The dipole moment functions 
converge against a unique function as the level of approx- 
imation is refined as shown in Fig. [3 indicating that the 
highest order CCSDT function is a good approximation 
to the physical dipole moment function. Furthermore, 
we have performed equivalent calculations on the isoelec- 
tronic molecules NaH and BeH+ 37., ,38^ ^39(1 to compare 
our results with other published calculations. The re- 
sults were in agreement within 5%, a level which is not 
critical for the simulations of the cooling schemes. The 
calculated Einstein coefficients are given in appendix 1X1 
We have now set up the model and acquired the pa- 
rameters entering the coupling matrix K in Eq. (Q and 
the solution can now be found numerically using stan- 
dard methods as described below. 



spectral distribution of the incoherent source used in the 
schemes of Fig.^b) andQIc) is found to be a square dis- 
tribution with the maximal density allowed on the rota- 
tional transitions from {v — 0, N = 1) ^^ {ly = 0, N = 2) 
to (ly ^ 0,N ^ 3) ^ {v = 0,N ^ 4). Furthermore, we 
showed that the spectral radiation density reaching the 
molecular ions from a realistic lamp is approximately 5 
times the spectral radiation energy density of BBR at 
300K. This has been included as a constraint in the op- 
timization of the incoherent field. 
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C. Solving the population dynamics 

We model the dynamics of the cooling on the test case 
of MgH+ by solving Eq. Q |40| ■ In the population vec- 
tor, P, of Eq. Q we use N„iax — 20 since the popu- 
lation of this and higher-lying levels is effectively zero 
during the cooling process. In addition, we omit the 
[v = 2, N) states in the cooling schemes if the second 
excited vibrational state is not coupled by laser fields. 
The radiation density W{uj) at resonance between lev- 
els not addressed by lasers has been calculated from a 
Boltzmann distribution at 300 K plus incoherent fields 
from lamps as described below. The pulsed lasers are in- 
cluded by saturating the pumped transitions, described 
in Sec. m at a repetition rate of 100 Hz. In the simula- 
tion this is done by redistributing the population in the 
involved ro- vibrational levels at the given repetition rate 
according to the degeneracy of the levels. All simulations 
are made with populations which are initially Boltzmann 
distributed at a temperature of 300K. The shape of the 
incoherent field is chosen so it maximizes the final popu- 
lation in the rovibrational ground state. 

All simulations are made with the most abundant iso- 
topes, in this case 24MgiH+ (79%). 



D. Efficiency of ^E cooling schemes 

In Ref . |2J| we found that the optimized radiation den- 
sity at intermediate timescales induces transitions up to 
and including the peak of the population distribution in 
BBR alone at 300K. Specifically, for MgH+ the optimized 



FIG. 3: Population in the ro-vibrational ground state of 
MgH^ (X^E) vs. cooling time for the Raman (a) and direct 
(b) schemes presented in Figs.^^b) and (c) with the optimized 
distribution of incoherent radiation (solid) , quartz- filtered dis- 
tribution (dashed) and no incoherent source (dash-dot). The 
same simulation using the scheme of Fig. H^a) is depicted for 
comparison (dot). 

While the cooling efficiency at a given time depends 
critically on the ability to filter out radiation addressing 
the heating transition in the low-frequency end of the 
distribution, it is only weakly depending on the sharp- 
ness of the filter in the high-frequency end. This is il- 
lustrated in Figs. El and 01 by the inclusion of a sim- 
ulation using a square incoherent field addressing the 
transitions (jy == 0, TV = 1) ^ (z/ = 0, TV = 2) to 
(i/ = 0, iV = 7) ^^ {ly — Q,N — 8) roughly correspond- 
ing to the cutoff frequency of a crystalline quartz win- 
dow |4ll. This distribution will be referred to as the 
"quartz-filtered" distribution below. Simulations of the 
evolution of ro-vibrational ground state population of 
MgH+ during cooling with the direct and Raman scheme 
using these two incoherent fields are presented in Fig. 
together with the results obtained without the inclusion 
of an incoherent source and those obtained by applying 
the scheme of Ref. |2l| (Fig.Uta)). 

For very short cooling times no significant difference 
between the schemes is seen as the relatively slow rota- 
tional transitions have not yet set in. On intermediate 
time-scales the effect of the added incoherent field is ev- 
ident and the o ptim ized scheme has an advantage to the 
scheme of Ref. ^l] at times less than ^ 100 s. At long 
times the slower depletion of the {v = 0, N = 1) state us- 



ing the incoherent field rather than a laser, as well as the 
heating effect of the added radiation makes the scheme 
of Ref. I23 more effective than the other schemes. 

The schemes presented here has the advantage of 
reaching significant cooling after ^ 30 s which, com- 
bined with the modest demands to coherent light sources, 
makes them experimentally attractive. Anticipated per- 
formance of traps for neutrals give storage times ex- 
ceeding 10 s, comparable to the timescale of the cooling 
schemes applied on MgH"*" Q and in the same regime as 
ArH+ , which is a faster candidate |23 ■ Hence the appli- 
cation of the schemes may be considered to create rota- 
tionally cold neutral molecules in the presence of BBR. 

The population distribution after 60 s of cooling is 
compared with the initial Boltzmann distribution in Fig. 
21 The depletion of the rotational levels above the N = 2 
"pump state" is evident. The difference between using 
the optimized and quartz-filtered spectral distribution of 
incoherent light can be seen in the figures, but the effect 
is very limited. 

The final population in the ro- vibrational ground state 
of just below 80%, c.f. Fig. 01 corresponds to the ground 
state population of a thermal ensemble of MgH+ at ~7 
K. 



TABLE I: Overview of quantum numbers describing the rovi- 
brational state of a molecule, neglecting nuclear spin. ^ de- 
notes a unit vector along the internuclear axis. 
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FIG. 4: Population distribution in MgH+ (X^E) after 60 s 
of cooling using the Raman (a) and direct (b) schemes with 
the optimized energy distribution of the incoherent source 
(black), the quartz-filtered energy distribution (gray), com- 
pared to the initial population distribution at 300 K (white) . 
The ground state population after cooling with the optimized 
incoherent field corresponds to that of a thermal ensemble at 
~7K. 



IV. COOLING SCHEME FOR MOLECULES 
WITH ROTATIONAL SUB-STRUCTURE 

In the previous sections, we discussed cooling schemes 
applicable to molecular ions with their ro- vibrational en- 
ergy levels determined by molecular rotation and vibra- 
tion only. This will be the case if the relevant electronic 
state has vanishing total spin and if the projection of the 



Label 


Definition 


L 


Total electronic orbital angular momentum 


A 


Projection of L on internuclear axis 


N 


Angular momentum of molecular rotation 


S 


Total electronic spin 


E 


Projection of S on internuclear axis 


Ms 


Projection of S on laboratory Z-axis 


n 


A + T, 


K 


Sum of AT and A • C 


J 


Total angular momentum of molecule neglecting 
nuclear spin 


Mj 


Projection of J on laboratory Z-axis 



orbital angular momentum of the electronic state along 
the internuclear axis is zero, i.e., in ^S states. We now 
turn to the other electronic ground states found in lighter 
diatomic hydride ions: ^E , '^E and ^11 

A range of quantum numbers will be needed to describe 
the rotational sub-states of the molecules to be discussed. 
We follow the notation of Herzberg |2J] designating the 
quantum numbers as indicated in Tabled The meaning 
of the coupled angular momenta is explained below. 

We now treat Hunds coupling case (a) and (b) sepa- 
rately and study cooling schemes for both cases. 



(2S+1) 



H-states; Hunds case (a) 



An interaction term of the form H'^° — AL ■ S will 
appear in the Hamiltonian if the projection on the in- 
ternuclear axis of both electronic spin, S, and electronic 
orbital angular momenta, L, arc nonzero. For moderate 
rotational excitations this will normally dominate over 
terms from the rotational Hamiltonian, H"^"^ = B ■ N^. 
It is therefore convenient to choose the Hunds case (a) 
basis set, consisting of basis functions In, S^^ J^MjAEJl) 
where n is collecting the quantum numbers defining the 
molecular state but not mentioned in Table ^ In this 
basis set, the unperturbed Hamiltonian, Hq is diagonal 
and the main perturbation term H^° is nearly diago- 
nal with the off-diagonal terms satisfying Afi = 0. The 
\n, S'^J'^MjAT.n) basis states are therefore a good ap- 
proximation to eigenfunctions with good quantum num- 
bers if |A| 3> -B ■ J- In the following section, we re- 
strict the calculation to the pure Hunds case (a) limit 
where this condition is fulfilled. •^H states are often close 
to this limit at low rotational excitations and they form 
the most interesting example of Hunds case (a) coupling 
for our purpose, as they are found as ground states of 
a number of molecules interesting for cooling, including 
NH+ and FH+ . 



1. Energy levels and selection rules 

The first order effect of H'^° is to split the electronic 
ground state into states according to the value of fi. For 
each of these states there will be a set of ro-vibrational 
sub-states arising from H^°* . 

In Hunds case (a), the molecule is well-described as a 
rotating symmetric top, for which the rotational energies 
are expressed by |24j 



F^{j) = B^{j{j + i)-n^). 



(9) 



Here J must take values greater than |rj — A^ | and the low- 
est rotational state will therefore, in general, have J 7^ 0. 
The overall structure of the molecular energy levels can 
be seen from the sketch of the modified cooling scheme 
inFig. Elfor 5 = \. 

The case (a) basis state in the laboratory frame can be 
written as a Wigner rotation of the corresponding wave 
function in the molecular rest frame 1431 
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{{r[}R\n)\ST.)Vij^^{al3^) 



(10) 



where {ri\,R ({r^},i?) are the electronic and internu- 
clear coordinates in the laboratory (body-fixed) frame. 
Finally Vf,j^{a/3'^) is an element of the Wigner rotation 
matrix evaluated at the given Euler angles, af3j AS] . The 
Honl-London factors S{J',J") are found as outlined in 
Appendix IBl 



S{J',J")^{2J" + l)x 

\{j"n"i{n'~n")\j'n')\^Ss'^s"S: 



(11) 
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where (J"n"l{n' - n")\J'n') is a Clebsch-Gordan coef- 
ficient. This result immediately gives us the following 
dipole selection rules 

AJ = 0,±1 but J^O^J^O 

AA = 0,±1, (12) 

AS' = AE = 0, 

which can also be combined to Ail — 0. 
2. Cooling schemes 



For ^n molecules we propose the cooling scheme de- 
picted in Fig. |31 where we distinguish between the pos- 
sible values of fi = ^ and ft — ^. Since only transi- 
tions with AJ = 0,±1 are allowed, we can pump pop- 
ulation from the first excited rotational state in the vi- 
brational ground state to the rotational ground state of 
the first excited vibrational level. The former is denoted 
the "pump state" in analogy with the nomenclature in 
Sec.m From the {i^ — 1,J = U) state spontaneous emis- 
sion brings population either back to the pump state or 



down to the ro-vibrational ground state. The cooling 
scheme must be applied for each populated f2 state indi- 
vidually. In Fig. we have assumed population of both 
fi = ^ and fl = |. In the absence of incoherent radi- 
ation this forms a pumping cycle where population ini- 
tially in the {v, J) — {Q,il + 1) state is transferred to 
the ro-vibrational ground state. As in the singlet case, 
the presence of BBR and possibly additional incoherent 
radiation from a lamp, will induce rotational transitions 
and thereby feed the pump state with population from 
higher-lying states. The entire population is therefore 
cooled. 

Cooling schemes for other Hunds case (a) molecules 
may be derived from straightforward generalization of 
the ^n scheme. 



J 



J 




II 
> 



O 

II 

> 




3/2 J- 



^n 



3/2 



FIG. 5: Cooling scheme for ll-states. Each of the two pos- 
sible values of Q results in a series of energy levels and it 
is necessary to cool the ^IIi and ^Ila separately if both are 
populated. Population is pumped from the first excited rota- 
tional state in the i' = vibrational ground state to the ro- 
vibrational ground state by a laser induced vibrational tran- 
sition and subsequent spontaneous decays. All the involved 
transitions are dipole allowed, cf., Eq. 11211 . Solid lines in- 
dicate laser pumped transitions while dashed lines indicate 
spontaneous decay paths. The A doubling is not shown in 
the figure. 



3. Numerical simulations 

The simulation is done using the approach described in 
Sec, nil CI but with the dipole transition matrix elements 
calculated using the Honl-London factors of Eq. Hll|l . 
We have chosen the molecule FH+ as an example of a 
^H ground state molecule. 



Since the spin-orbit coupling parameter A 
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cm~ is much larger in magnitude than the rotational 
constant 5 = 17 cm~^ FH+ is best described in the 
Hunds case (a) scheme |44^. The appropriate cooling 
scheme is depicted in Fig. |S| although it should be noted 
that, for FH+ , fi = | is the lower state. To model the 
cooling scheme we use the dipole moment functions in 
Rcf . 45] and the accurate spectroscopic data of Ref. [ij . 

In the cooling scheme of Fig. [S] the pumping is done 
from the first excited rotational level. This fact, com- 
bined with a large permanent dipole moment and hence 
rotational transition rate of FH"*" (2.57 Debye), makes the 
effect of the broadband incoherent radiation marginal. 
We have therefore performed the simulations without the 
inclusion of an incoherent source. The results of simula- 



tions are given for both the ^11 
Figs. El and 



1/2 



and ^n 



3/2 



states in 



Further splitting of the levels indicated in Fig. [S] will 
appear due to A doubling. The effect is largest in the 
^n 1 state where it has a magnitude on the order of 10 
GHz, which is more than one can expect to cover with the 
bandwidth of a single pulsed laser. Therefore the laser 
transitions indicated for the ^11 i scheme needs to be 

2 

divided into two. The splitting of the lowest ^H a state 
is an order of magnitude smaller, so it is not necessary to 
split that laser transition if a pulsed laser system is used. 
This leaves us with three laser frequencies to use for the 
cooling scheme if we assume that both ^ — ■^ and f^ = § 
are populated. 
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FIG. 6: Cooling efficiency of FH+ (X^H) in the two Q. sub- 
states. In the dipole approximation the two sub-states are 
uncoupled if pure Hunds case (a) applies. The cooling scheme 
will therefore be significantly simplified if one can design the 
experiment such that only the lowest f^ = § state is populated 
in the cooling scheme 



Complications arise if we are not in the pure Hunds 
case (a) scheme. This occurs if the rotational part of the 
Hamiltonian cannot be neglected compared to the spin- 
orbit part. Treated in the case (a) basis, the rotational 
part will produce non diagonal perturbations |43l- This 
would allow a coupling from {v = 1, J = fi) — » (v — 
0, J = O -I- 2) (the introduction of quadrupole couplings 
would have a similar effect). We do not expect this effect 
to be significant given the difference between |^| and B. 
We did, however, check the stability of the scheme when 
introducing such couplings and found that due to the fast 
rotational redistribution rates, the population that was 
coupled out of the cooling cycle by A J = 2 transitions 
would rapidly be taken back. The negative effect of such 
couplings is small (less than 10% decrease in cooling ef- 
ficiency) if the A J = ±2 couplings are less than 20% of 
the A J = couphng strength. 

It should be noted that since the coupling between the 
n states is also absent in the pure case (a) coupling, it 
would be possible to prepare the sample so that only the 



n 



sub-state is populated, due to its significantly 



lower energy. This would make the lasers addressing the 
other level superfluous. In that case only a single laser 
frequency is needed to cool the molecules. 
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FIG. 7: Cooling efficiency as function of cooling time for the 
two Q sub-states of the ^11 electronic ground state of FH^ . 
The cooling is seen to reach steady-state after < 10 s without 
the inclusion of an incoherent source, largely due to the large 
permanent dipole moment of FH^ . 



B. '^^•^+^)E-states; Hunds case (b) 

li B > \A\ or at high rotational excitations the Hunds 
case (a) basis functions will no longer be approximate 
energy eigenfunctions. If H^°* dominates, the Hunds 
case (b) basis, |nJ^MjA/"^S^A) is convenient as the total 
Hamiltonian is nearly diagonal in this basis. In particu- 
lar this is fulfilled for ■^'^+^S states which are common as 
electronic ground states of light diatomic molecular ions, 
including BH+ (X^E) and 0H+ (X^S). Below we treat 
the ^S and "^S cases separately. 



1. Energy levels of doublet states 

The sub-states of a rotational level in a molecule in a 
^E state are split due to the interaction of the spin of the 
unpaired electron and the molecular rotational angular 
momentum. This is due to the spin-rotation Hamilto- 
nian, H"^ = jN ■ S, with 7 denoting the spin-rotation 
coupling constant. The resulting energies of the doublet 
are given by |2J| 

FiiN) = BN{N + l) + ^jN, (13) 

F2iN) - BNiN + l)-^-/{N + l), (14) 

and the sub-states are denoted Fi and F2 ior J = N + ^ 
and J = N — ^ respectively. 




FIG. 8: Cooling scheme for ^E-states. Due to the spin- 
rotation coupling each rotational quantum state A'' split into 
two sub-levels with J=|iV + i|,|iV-i|. The dipole-allowed 
vibrational transitions are indicated on the figure using solid 
lines for laser pumped transitions and dashed lines for subse- 
quent spontaneous decay paths. 



2. Energy levels of triplet states 

Molecular ions in ■^E electronic states will, apart from 
the spin-rotation splitting discussed above, have an ad- 
ditional splitting from the coupling of the electronic spin 
of the two unpaired electrons. Such states are relatively 
rare, as pairing of the electronic spins is usually favored. 
Nevertheless, the ionic hydrides in the I6th group of the 
periodic table, including OH"*" and SH"*" , have such elec- 
tronic ground states and we therefore consider the appli- 
cability of the cooling schemes to such states here. The 
energies of the three spin sub-states are given by |24| 



3. Selection rules 



Fi{N) 
F2iN) 
F3{N) 



BN{N 
BN{N 
BN{N 



1) 



1) 



2A(A^ + 1) 


2iV + 3 


2AA^ 



2N - 1 



+ 7(7V + 1),(15) 
(16) 

jN. (17) 



In analogy with the doublet case Fi, F2, and F^ denotes 
the sub-states with J ^ N + I, J ^ N and J ^ N - 1, 
respectively. In the expression 7 is the spin-rotation cou- 
pling constant and A is the spin-spin-splitting constant. 
The latter is normally an order of magnitude or more 
larger than 7 and, consequently, the multiplet splitting 
of triplet states at moderate rotational excitations are 
much greater than the corresponding splittings of a dou- 
blet electronic states. 



N 




FIG. 9: Cooling scheme for E-states. Due to spin-spin and 
spin-rotation coupling each rotational quantum state N split 
into sub-levels with J = \N + 1\,N, \N - 1|. The dipole- 
allowed vibrational transitions are included in the figure with 
solid lines to indicate laser pumped transitions and dashed 
lines to indicate spontaneous decay paths. 
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In Hunds case (b) the good quantum numbers are 
N, S, J, Mj and A. We therefore write the eigenfunctions 
in the laboratory frame as 



{{r,},R\nJMjNS,K) = ^ ^^^x 

S N 

J2 Y. {{rl},R\n){NMNSMs\JMj) 



(18) 



-N 



SMs)V%,.{af3^l 



where {ri},R {{r[},R) are the electronic and internu- 
clear coordinates in the laboratory (body-fixed) frame. 
We then follow the approach of Appendix El to get the 
Honl-London factors 



S{J', J") = {2N" + 1)(2J' + 1)(2J" + l)x 

.2 {S' N" J"^' 

1 J' N' 



(N"A"1{A' ~A")\N'A'y 



Ss'S", 



(19) 



S' N" J" 

where < ^ ^ > is a 6j symbol 43}. The following 

selection rules are extracted 
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FIG. 10: Population in the lowest rotational state of 
BH^ (X^E) as function of cooling time. Simulations are made 
using the scheme of Fig. |H| with BBR only (dashed line) and 
with the inclusion of the field from an incoherent source ad- 
dressing the A'^ = 1 ^ TV = 2 and A = 2 ^ A = 3 transitions 
(solid line). We see that a significant improvement is obtain- 
able using the incoherent source. In line with our experience 
from MgH"*" there is only a couple of percent loss of cooling 
efficiency when using a softer low-frequency pass filter, for ex- 
ample, letting the broadband incoherent radiation extend to 
include transitions up to and including A = 4 ^ A = 5. 



AJ = 0,±1 but J = O^J = 

AN = 0, ±1 but AiV 7^ if A' = A" 

AA:=0,±1, 

AS-^O. 



(20) 



4- cooling scheme 

The cooling scheme proposed for Hunds case (b) 
molecules closely resembles the singlet cooling scheme. 
It is depicted in Fig. |S1 for ^S molecules and in Fig. |5| 
for ■^E molecules. The optical pumping is done from the 
(i/ = 0, A^ = 2, J) set of states to {ly = 1,N = I, J'). 
Then dipole allowed spontaneous decay will result in 
transitions back to the "pump states" or to the non- 
degenerate ro- vibrational ground state. The only change 
to the scheme when compared to the singlet case is to as- 
sure the addressing of all sub-states in the A^ multiplet. 
This is possible because the AA^ = ±1 selection rule for 
S states from Eq. (I2()|l is the same as in the singlet case. 
The role of BBR and additional incoherent radiation is 
the same as in the previous schemes. 

The number of transitions to be pumped is three for 
the ^S states and two for the ^E states. The splitting of 
the levels in the former is expected to be much larger than 
for the ^E case, since the spin-spin coupling parameter 
A is much greater than the spin-rotation parameter 7 as 
mentioned in Sec. lIVBl 




12 3 4 5 6 7 
Rotational level 

FIG. 11: Population in the lowest rotational states of 
BH^ (X^E) after cooling in 120 s using the incoherent ra- 
diation from a lamp addressing the A = 1 -^ A = 2 and 
A = 2 ^ A = 3 transition (Black columns) and in BBR 
only (grey columns). The initial 300K Boltzmann distribu- 
tion is included for comparison (unfilled columns). We note 
that slightly better cooling efficiency should can be obtained 
using longer cooling times, cf Fig. IIUI This is, however, im- 
practical and the obtainable improvements would be rather 
small. The substructure of the rotational levels is included in 
the simulation but omitted on the figure. 
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5. Numerical simulations; BH^ (^'^ ) 

Here we treat BH+ as an example of a ^S ground state 
molecule and discuss the molecule specific parameters 
and their implications on the cooling schemes. 

The numerical simulation is done for ^^B^H+ which is 
the dominant isotope (80%). We use the potential en- 
ergy and dipole moment functions of Ref. [43 ■ With 
those functions, we use the approach of Sec. IIV B 31 to 
calculate the matrix of Einstein coefficients between ro- 
tational and vibrational states. Finally, we make the sim- 
ulation as described in Sec. IlllCl but with the modified 
energy level structure. If one neglects fine-structure, the 
laser wavelength for the two, then identical, pump tran- 
sitions depicted in Fig. |S1 is f^o = 4.17^m. The real 
resonant transition frequencies are shifted from this cen- 
tral frequency through Eq. \i'A\ where 7 = —0.014 cm^^ 
|24|. This gives a splitting of laser frequencies, including 
fine structure, of ±0.007 cm"! ~ ±210 MHz. This differ- 
ence is comfortably smaller than the typical bandwidth 
of a pulsed laser system. The hyperfine coupling coeffi- 
cient has, to our knowledge, not been calculated. Typi- 
cal values are, however, on the order tens to hundreds of 
MHz, allowing us to address all hyperfine substates with 
the same pulsed laser system. Hence, it is reasonable to 
expect that for practical implementations only a single, 
pulsed laser frequency is needed. 

The results of a numerical simulation are given in 
Figs. 1101 and 111! We note that the convergence is quite 
slow compared to what we saw from MgH"*" and FH+ . 
Optimal cooling is not obtained until after ~ 2 min- 
utes. This is not too critical as 60 % of the popu- 
lation is in the ground state after 20 s. As expected 
from the discussion in Ref. [22| we find, that the op- 
timized distribution of the incoherent source addresses 
the transitions (i/ = 0, TV = 1) ^ (i/ = 0, iV = 2) and 
(jy = 0,iV = 2) ^ (1/ = 0,iV = 3). Similarly it is con- 
firmed, that the cooling efficiency has little sensitivity 
towards the high frequency cutoff of the incoherent field. 



6. Num.erical simulations; OH^ f'^) 

As an example of a molecule with the "^E ground state 
we have chosen 0H+ . This molecule plays an important 
role in the chemistry in comet tails 1 131 . the upper earth 
atmosphere and interstellar clouds |43|. The electronic 
ground state of 0H+ is '^E^ . The effect of hyperfine split- 
tings is expected to be much smaller than the bandwidth 
of a typical pulsed laser system due to the nuclear spins 
/ = and / = i of O and H respectively. Hence the 
molecule is well-described by the level scheme of Fig. 
The frequencies of the three laser beams required are 
found from Eq. IjlSII and the constants 7 — —0.0147 cm~^ 
and A = 2.13 cm~^ uM- The wavelength of the un-split 
transition is 3.3 /xm with the three sub-transitions shifted 
-3.2 GHz, and -12 GHz with respect to it. This splitting 
is to large to be covered by a single broad laser unless one 
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FIG. 12: Top: Born-Oppenheimer potential curves for X E 
OH^ calculated by Gaussian with various theoretical models 
compared to the calculation of Ref. 1431 . All our calculations 
are done in a 6-311-|--|-G basis set except the solid black line 
which is made in the generally more accurate aug-cc-pVTZ 
basis |31|. The curves agree close to the equilibrium, 1.03A, 
for the MP4 (Fourth order M0ller-Plesset perturbation the- 
ory) and Coupled Cluster approaches (CCD, CCSDT) indi- 
cating an accurate level of theory. The different methods are 
described in Refs. HHH HI ES IH ■ Bottom: Dipole mo- 
ment function calculated with Gaussian using similar levels 
of theory and basis sets. The agreement between the calcula- 
tions is reasonable and the effect of using the larger basis set 
for the CCSDT theory is not visible on the given scale, but 
our results show some discrepancy with the results of Ref. |4^ . 
This small discrepancy, however, has very little effect on the 
cooling scheme. The classical turning points for the vibra- 
tional ground state are marked on the common abscissa at 
0.95 and 1.15 A. The dipole moment functions are given in 
center of mass coordinates. 



finds a way to generate shorter and thereby broader and 
more intense pulses in this wavelength regime. This is an 
obvious experimental complication that will often arise in 
the case of '^E states due to the generally large value of 
the spin-spin splitting constant A. It should, however, be 
noted that the 3 GHz may be covered by a single pulsed 
laser, leaving only two laser frequencies in the cooling 
scheme. We have calculated the dipole moment func- 
tions of OH"*" and compared our results to Ref. [ig in 
Fig. El In the simulations we use the function obtained 
in the CCSDT (aug-cc-pVTZ) calculation. 
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The final population distribution in the numerical sim- 
ulation is given in Fig.^| The scheme is both faster and 
more effective than what was found for MgH+ . This 
can be understood from comparison of the Figs. |21 and 
[T^ A larger gradient of the dipole moment function of 
0H+ results in a higher effective pump rate from N = 2 
to N = 0. As with MgH+we see a significant increase 
in the cooling efficiency when introducing broadband ra- 
diation from an incoherent source to deplete the A^ = 1 
population. 

The simulation shows the efficiency of the rotational re- 
distribution in the ^S state. Considering the nonzero line 
strengths for transitions between the Fi,Fj, (i ^ j) se- 
ries of states, provided AA^ = ±1, one could be tempted 
to omit one or more laser frequencies expecting rota- 
tional redistribution to empty the remaining substates 
by rotational transitions through neighboring A^-levels. 
Unfortunately such redistribution rates, requiring two 
or more rotational transitions through specific substates, 
are much too slow to have a significant effect on the cool- 
ing scheme. Therefore each of the three laser frequencies 
are needed to make the cooling scheme effective. In ac- 
cordance with the previous results we find that the op- 
timized distribution of incoherent radiation from a lamp 
addresses only the {ly ^ 0,N = 1) ^ {f = 0,N = 2) 
transition. 

Finally, it should be noted that ^■^+^I]-states are al- 
ways cases of pure case (b) coupling due to the vanishing 
orbital angular momentum and the selection rule in A^ is 
close to exact. This stands in contrast to "^^+^11 states 
which often have effects of intermediate coupling which 
will complicate the suggested case (a) cooling scheme fur- 
ther. 
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Rotational level 



FIG. 13: Population in the lowest rotational states of 
OH"*" (X^E) after cooling in 10 s using the incoherent radi- 
ation from a lamp addressing the N — 1 ~* N — 2 and 
A = 2 — > A = 3 transition (Black columns) and in BBR only 
(grey columns). The initial 300K Boltzmann distribution is 
included for comparison (unfilled columns). The spin sub- 
structure of the rotational levels is included in the simulation 
but omitted on the figure. 



V. SUMMARY 



We have presented cooling schemes for rotational cool- 
ing of translational cold molecular ions in the ^E , ^E , 
■^E and ^11 electronic ground states. For all but the rela- 
tively rare ^E electronic state the schemes can be realized 
by optical pumping with a single pulsed laser beam, pos- 
sibly combined with the inclusion of a broad-band inco- 
herent source. They are therefore experimentally attrac- 
tive, and preliminary experiments are presently under 
way with MgII+ . 

Possible applications include high-precision spec- 
troscopy and measurements of absolute reaction rates 
with molecular ions in a single, well defined quantum 
state. This could for example be used to study disso- 
ciative recombination with unprecedented resolution or 
molecular reactions in interstellar media or comet tails 
[ill Il2 |. Ultimately, the access to cold molecular ions 
could be used in implementations of quantum logics. 
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APPENDIX A: EINSTEIN COEFFICIENTS 
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TABLE II: Einstein coefficients for selected transitions in s ^ and corresponding transition frequencies in cm ^. In the table 
the quantities have the following meaning: ^4''°* = A{u = 0, N = 1 ^ u = 0, N ^ 0), A''''' = A{u = 1, N = 1 ^ u ^ 0, N ^ 0) 
for the E-states and A''"* = A(iy = 0, J = ft + 1 -^ u = 0, J = ft), A"''' = A{iy = 1, J ^ n ^ ly = 0, J = fl) ior the H-state. A 
similar notation is used for the transition frequencies. The pure rotational transition rates (first column) indicate the rotational 
redistribution speed while the vibrational transition rates give the spontaneous decay rate from the excited vibrational state 
in the pumping schemes. The data largely explains the qualitative difference in cooling efficiency for the molecular ions. Large 
rotational redistribution rates indicate a fast scheme, while fast spontaneous decays from the excited vibrational state indicate 
high effective pump rate, i.e. high cooling efficiency. The data for BH+ is found using the data from Ref. 47] and the computer 
programme of Ref. |30| from which the Honl-London factors were corrected to conform with the multiplet expressions of Sec. llVI 
Data for OH+ was found using the same approach and data from Fig. 1121 Finally the data on FH+ was obtained using the data 
of Refs. Ellii. 
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APPENDIX B: HONL-LONDON FACTORS 

For completeness we include the details of the deriva- 
tion of Eqs. iIHIIjCU and itT^ . 



Which is the same as Eq. (0J except for the degeneracy 
factor. We now move to a molecule-fixed coordinate sys- 
tem. We define the electronic dipole moment function 
by integrating the dipole operator, 7Vi'™°', over the elec- 
tronic variables, Te 



1. 'E ground state 



Eq. I0J must be modified if 'i'm and ^„ are degen- 
erate. The effective Einstein B-coefficient is found as 
Bn,m = X)u Sf "''''"^ , where /i, £, denote the sub-states 
of ^„ and ^„i , respectively, and gn the degeneracy of the 
initial (upper) state. One summation is done to include 
transitions to all sub-states of the final state, m, while 
the remaining terms correspond to averaging the result 
over the sub-states of the initial state. We then define 
the total transition dipole moment for degenerate states 
as 



|A^,: 



E 



^n^ ,r. 



(Bl) 



where the summation is done over all transitions between 
sub-states of the system. The Einstein coefficients be- 
tween degenerate states then take the form 



A. 



AT^n,r. 






(B2) 



:Br, 



Dr\R)= M{rl},RrM"^'''M{K},R)dT,- (B3) 



Here we stay in the electronic state defined by the wave 
function tpe{R) in the body-fixed frame. We have calcu- 
lated D^°'-{R) ah initio with Gaussian [23. Details of 
these calculations are molecule-specific and will be given 
below. 

To transform the dipole moment to the laboratory 
frame we now specialize to ^S states, postponing the 
general solution to Sec. II VI For ^S diatomic molecules 
the cylindrical symmetry of the potential will ensure that 
Z>™°'(i?) points along the internuclear axis. Hence the Z- 
component of D'^°^{R) in the laboratory system is given 
by 



D'^'"'{R)z = r>™°'(i?) ■ Z = D'^°'{R)cose 



(B4) 






The molecular states are degenerate so it is necessary 
to sum over all sub-states to obtain the transition dipole 
moment defined in Eq. (|Bip . In the case of ^S molecules 
this corresponds to summing over all projections Mj of 
the molecular angular momentum J . In carrying out the 
summation over the sub-states in Eq. (|Bip the selection 
rule AMj = 0, ±1 makes it possible to rewrite the expres- 
sion as a single sum over /.t = Mj which can be related to 
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the total transition dipole moment. Since the transition 
probabihty must be independent of the orientation of the 
laboratory coordinate system we have 



Mj 



b\2 



(B5) 



Mj 



Inserting Eq. IJB4(1 in Eq. ^ we find 



j-\lab / ^i,mol ( Ci A. T^\* j^Tfiol { 



V'"%„ (61, 0, RYD'^°\R) cos 61 X 
V'r"^ (^, 0> R)R^ sin edRdOdct), 



(B6) 



where ipr„,iyn and iprm,^^ are the remaining ro- vibrational 
wave functions obtained after the integration over elec- 
tronic coordinates in Eq. I|B3I) . Now, we assume that the 
ro-vibrational wave function may be written as a prod- 
uct, ^r-„^^JR,e,<j)) = ^rSQA)UAR)- Then 



\D 



with 



lab\2 



L 






UjR)Dr\R)U^iR)R'dR 



(B7) 



L 



JtuMtt 



^lab 



$r {e,<j))cose<i>';"'{e,(t))sm9ded(t). (bs) 



Combining the above equations and performing the in- 
tegral over Euler angles, while writing the Wigner ro- 
tation functions as an expansion over Clebsch-Gordan 
coefficients '51^, one finds the dipole moment transition 
matrix elements (fc — 1) 



{n'J'M'j\Tli{r,},R)\n"J"M'J)^ 

{j"M'jii\j'M'j) {j"n"kA\j'n') . 



(B13) 



Summing over the projections of J and emission direc- 
tions one finds the line strength 



J2 I {n'J'M'jl Tl{{n}R) \n"J"M'}) p = 

M'j.M'J 

{2J" + l)\{n'v'\Tl\n"v")\^x 

{j"n"i{n' - Q.")\j'n') \^6{s', s")5{y,\ s"). 



(B14) 



Finally we find the Honl-London factors in Hunds case 
(a) 



5(J',J")-(2J" + l)x 



(B15) 



Defining 



S. 



J,^,J,^ 



M„ 






(B9) 



known as the Honl-London factors 24, 27, '50]. We com- 
bine the above results with Eq. (|B5(I to find the total 
transition dipole moment entering Eq. (|B2p 



|X>, 



Sj„^j„ 



U„iR)D,{R)UjR)R^dR 



(BIO) 



2. Hunds case (a) 



We use the Hunds case (a) eigenfunctions in the lab 
frame from [H (cf. Eq. (HUJl) 

{{r.,}R\nJMjnS^) = 

2/11 . (Bll) 

^^({ra,i?|n)|5E)P]J,^^(a/?7), 

and write the ^th component of the A;th moment transi- 
tion operator in the laboratory frame, T^'"', as a similar 
rotation of the operator working in the molecular rest 
frame 



3. Hunds case (b) 

We gave the Hunds case (b) eigenfunctions in the lab- 
oratory frame in Eq. (|18|) 



{{n}R,nJMjNMNSMs) = y ^^^J^x 

S N 

E E {{rl},R\n){NM^SMs\JMj)x ^^^^> 



Ms = -S Mn = -N 



\SMs)VZj^^{aM)- 

The rotated dipole moment operator was given in a gen- 
eral form in Eq. ljB12p . We then use the identities |43 



T^l'Dl,,= E {nmkl\N'M')x 

N'M'f,' (BIT) 



and 



T^imPx^i.'^^ — ^TTVT'"'^-^"™'^ 



6L\6ra.,^lh,K., (B18) 



f-2ir , f-2TT , i-TT 



where J dil = L da L d'y L dj3sin[3. One thereby 



Tl'{{ri}R) = E T'l{{r[},R)'D'l'^{a(3^). (B12) finds the expression for the dipole matrix element 



A=-fc 



15 



{n'J'M'j\Tl{{r[},R)\n"J"M[i) 



2N" 



2N' 



{n'u'\Tl,_^„ \n"u") {N"k"l{K' - k")\N' k') x 






{N"M';,ll\N'M'j,) (N'M'^S'M'slJ'M'j) {N" M'^^S" M'^\J" M'}) Ss'S"Sm'm'4 



(B19) 



This is summed over the projections of J and squared 
to find the dipole transition probabihty. The task is sim- 
pUfied by rewriting the products of Clebsch-Gordan co- 
efficients in terms of Wigner 6j symbols [53 • After some 
algebra one then finds 



{n'J'N'\Tl\n"J"N")\^ 



1 



{2N" + l)x 



(2J' + 1)(2J" + 1) (nV'l Tl,_j^„ \n"v"f x 
(iV"A"l(A'-A")|iV'Ay|^ ^'' ^'1 . 



(B20) 



Hunds case (b) 



5(J',J") = (2iV" + l)(2J' + l)x 
(2J" + 1) (Ar"A"l(A' - A")\N'K'f x 

is N" J"l' 
1 J' TV' ^"'^"• 



(B21) 



Summing over the emission directions cancels the factor 
of o, leaving the expression fo r the Hnl-London factor in 
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